Pion photoproduction in a dynamical coupled-channels model 
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The charged and neutral pion photoproduction reactions are investigated in a dynamical coupled- 
channels approach based on the formulation of Haberzettl, Huang, and Nakayama [Phys. Rev. C 
83, 065502 (2011)]. The hadronic final-state interaction is provided by the Jiilich nN model, which 
includes the channels nN and r]N comprising stable hadrons as well as the effective ti-kN channels 
ttA, gN , and pN . This hadronic model has been quite successful in describing nN — s- -kN scattering 
for center-of-mass energies up to 1.9 GeV. By construction, the full pion photoproduction current 
satisfies the generalized Ward-Takahashi identity and thus is gauge invariant as a matter of course. 
The calculated differential cross sections and photon spin asymmetries up to 1.65 GeV center-of- 
mass energy for the reactions 7p — >■ n'^n, 7p — > tt^p, 771 — >■ 7r~p and 771 — >■ yr'^n are in good agreement 
with the experimental data. 
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I. INTRODUCTION 

Presently, there is intense experimental effort to study 
the production and the decay of baryon resonances from 
threshold to invariant collision energies of about = 2.8 
GeV, to obtain information about the non-perturbative 
sector of Quantum Chromodynamics, as discussed, e.g., 
in the recent review To deduce resonance param- 
eters from the experimental data, one has traditionally 
relied on partial- wave analysis. Most of our knowledge on 
baryon resonance masses is due to the classic partial- wave 
analyses of pion-nucleon scattering by Cutkosky 0, 
Hohler and Arndt and their collaborators. In 

the energy range under consideration, however, reaction 
channels other than pion-nucleon open up, which leads 
to ambiguities in the analyses. Moreover, the number of 
partial waves required scales with the energy. This sit- 
uation calls for other theoretical approaches. There are 
a few physical principles that an analysis should respect, 
such as unitarity and analyticity of the S-matrix, and, 
as a matter of course, gauge invariance if photoproduc- 
tion is considered. The various theoretical approaches 
to coupled-channels problems can be grouped into three 
broad classes, (i) unitarized chiral perturbation theory, 
(ii) dynamical coupled channel approaches, and (iii) K- 
matrix approaches. All classes guarantee the two-body 
unitarity of the 5-matrix by deriving the T-matrix from 
a Bethe-Salpeter or Lippmann-Schwinger equation, for- 



mally T — V + VGqT, where V denotes the scattering 
kernel defined by a set of Born diagrams, while Go stands 
for the intermediate two-particle propagator. The three 
classes differ by the choice of the scattering kernel V and 
the propagator Go- 

(i) In photoproduction reactions, the threshold re- 
gion is understood in terms of chiral perturbation the- 
ory (ChPT) 0-[ll| (for a recent comprehensive review 
see [ill). Unitarization of the interaction allows one to 
extend the applicability into the resonance region. Such 
chiral unitary approaches respect chiral symmetry, a fun- 
damental property of the strong interaction. Moreover, 
the scattering kernel V is obtained within a systematic 
counting scheme which limits the number of admissible 
diagrams according to the order of the chiral expansion. 
A scheme to gauge-invariantly couple the photon to the 
unitarized amplitude was developed in Ref. Il3ll and ap- 
plied to kaon and eta photoproduction [i3, llSl ]. See also 
Ref. [3l for a gauge-invariant unitary framework in the 
context of the recently discovered structure in photo- 
production on the neutron. For earlier works on pho- 
toprocesses in the chiral unitary framework, see, e.g., 
Refs. fTMi. 
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Chiral unitary approaches usually concentrate on 5*- 
waves although some authors consider higher-order terms 
in the scattering kernel, thus allowing them to study the 
S and P partial waves within this approach [2l| . 
A unitary coupled-channels model was developed in 
Ref. [13 where the partial-wave amplitudes for the 
and TriV states are obtained by analytic extrapolations 
of the subthreshold reaction amplitudes computed in 
ChPT. Chiral unitary approaches respect the analytic- 
ity of the T-matrix by keeping both the real and the 
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imaginary parts of the two-particle propagator Gq. An 
interesting physical consequence of analyticity is the pos- 
sibility to generate bound states or resonances by meson- 
baryon dynamics alone. The A^*(1535), A*(1700), and 
other resonances have been claimed to be dynamically 
generated psI - Esj . While chiral unitary approaches are 
certainly the most elegant ones theoretically, their actual 
applications to the analysis of data and the study of res- 
onances has been limited to low partial waves and an 
energy range well below 2 GeV. 

(ii) Dynamical coupled-channels approaches employ ef- 
fective meson-baryon Lagrangians to define the scatter- 
ing kernel V of the Bethe-Salpeter equation, giving up the 
chiral counting scheme for V. Analyticity of the T-matrix 
is guaranteed by solving the Lippmann-Schwinger equa- 
tion employing the full two-particle propagator. There 
is no restriction on the partial waves, which is essential 
for analyzing observables at higher energies. The scat- 
tering kernel V includes the exchange of mesons in the t- 
channel and the exchange of baryons in the w-channel and 
therefore generates correlations between different par- 
tial waves. As such approaches are a Lagrangian based, 
SU(3) symmetry allows also to correlate different reac- 
tion channels. Those correlations introduce an energy 
dependence to the non-resonant background which is not 
due to the presence of s-channel resonances. In principle, 
the method can generate baryon resonances dynamically, 
a feature shared with the chiral unitary approach, but 
in actual calculations, in most cases resonant s-channel 
driving terms are required for quantitative rendering of 
the data because the partial-wave correlations put tight 
constraints on the free parameters of the model that in- 
hibits to some degree the possibility of generating reso- 
nances dynamically. Still, it should be stressed that due 
to the non-resonant background, only a small number of 
s-channel resonances is required which should prevent the 
approach from claiming spurious resonances in the data 
[29|. As dynamical coupled-channels approaches respect 
analyticity, the analytic continuation of the amplitudes 
is possible and provides the poles and the residues of 
the S'-matrix on the various Riemann sheets. As field- 
theoretical quantities, the poles of the ^-matrix do not 
suffer from ambiguities, as for example Breit-Wigner fits 
do, and thus provide a more model-insensitive way of 
characterizing the baryon resonances. At this point, we 
would like to stress the importance of a dynamical treat- 
ment of three-body cuts, as they appear, e.g., in interme- 
diate ttttN states. In Ref. [1^ it has been shown that the 
resulting branch points in the complex plane are impor- 
tant for the extraction of the baryon spectrum; in models 
that do not contain these points, they may be simulated 
by poles, leading to potentially erroneous results. 

An actively developed dynamical coupled-channels ap- 
proach is the Jiilich model [sol-ls^l that includes the nN, 
r}N , aN, pN, and ttA hadronic channels. This ap- 
proach provides the hadronic part of the interaction in 
the present study, as specified below. Quite recently, the 
Jiilich model has incorporated the KA and KY, chan- 



nels '345, using SU(3) symmetry, and a global fit to the 
corresponding experimental data is in progress. A fea- 
ture in the Jiilich model, unique among modern dynami- 
cal coupled-channels approaches, is the treatment of the 
important t-channel exchanges with p and a quantum 
numbers: those are largely fixed using crossing-symmetry 
and dispersive techniques from NN — > tttt data [30| . 
thus drastically reducing the model dependence at this 
point. Progress in a different direction has been recently 
achieved in Ref. |33] where it has been shown for the 
Jiilich and similar models that the approach is suited to 
analyze upcoming lattice data on the baryon spectrum. 
A scheme has been developed to address finite-volume 
effects and lattice levels could be predicted. 

Among the various groups developing dynamical 
coupled-channels approaches [H, [s^l is the EBAC group 
[38| . Their approach includes explicit nnN contribu- 
tions in an approximate manner in an effort to satisfy 
some aspects of three-body unitarity. Like most coupled- 
channels models, it includes the channel in the one- 
photon approximation [395 . It is not gauge invariant, 
however. The EBAC approach has been applied exten- 
sively [39I - I4H in the analysis of pion photo- and electro- 
production data. 

(iii) Another class of approaches used to analyze pion- 
and photon-induced reactions are if -matrix models. Dy- 
namical coupled-channels approaches require a price for 
what is delivered: one has to solve coupled integral equa- 
tions. The technical effort can be reduced by approxi- 
mating the two-body propagator Gq in the Lippmann- 
Schwinger equation: omitting the real part of the two- 
body propagator, one reduces the integral equations to a 
set of algebraic equations. In such if-matrix approaches, 
unitarity is still respected due to the presence of the imag- 
inary part of Go, but analyticity is lost. As a conse- 
quence of the approximation, only on-shell intermediate 
states are taken into account when solving the scattering 
equation, while the principal-value (dispersive) parts of 
the scattering equation are neglected, which suppresses 
the full contributions of the virtual two-body interme- 
diate states and in general results in a reduction of the 
strengths of multiple-scattering contributions. The res- 
onance parameters may compensate for this approxima- 
tion. Due to its technical simplicity and flexibility, the K- 
matrix approach has made possible the quantitative re- 
production of a large body of experimental data. Recent 
applications of the method can be found in Refs. [7„,8i,,42i- 
|46| . A variation of the standard i^T-matrix approach (uni- 
tary isobar model) was developed by the Mainz group 
47, 48]. 

In the present work, we present a gauge- invariant treat- 
ment of pion photoproduction. A feasibility study of this 
reaction employing the Jiilich TriV model was presented 
in Ref. [1^ . The corresponding differential cross sections 
for neutral and charged pion photoproduction reactions 
were found to be in reasonable agreement with the data 
up to the total center-of-mass (cm.) energy oiW = 1.25 
GeV. In this study, we will extend and refine the calcu- 
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lations of Ref. |49| employing a more recent version of 
the Jiihch model [3j| for the hadronic part of the am- 
plitude, and combine it with a novel way of solving for 
the photoproduction current in a gauge- invariant manner 

The Jiilich ttN model [3TI - [33t is based on time-ordered 
perturbation theory (TOPT) fsi']. It is a coupled- 
channels meson-exchange model including the ttN and 
77 TV channels as well as the ttA, aN, and pN effec- 
tive channels which implicitly account for the resonant 
part of the ttttN channel. The interaction kernel corre- 
sponding to the t- and u-channel diagrams is constructed 
based on the (effective) chiral Lagrangians of Wess and 
Zumino [s^ . [53| , taking into account the corresponding 
lowest-order non-vanishing terms. This Lagrangian is 
supplemented by additional terms for the coupling of 
A, w, rj, ao, and a [s^ [s^. For the s-channel dia- 
grams, apart from the bare nucleon pole contribution 
which is renormalized by the coupling to the nN contin- 
uum state to reproduce the physical nucleon, the inter- 
action kernel includes eight genuine resonances, namely 
5ii(1535), 5ii(1650), 53i(1620), P3i(1910), Pi3(172oj, 
£)i3(1520), P33(1232), and 1)33(1700). The bare gen- 
uine resonances get their dressed masses and widths from 
the re-scattering of baryon-meson continuum states; the 
Pii(1440) (Roper) resonance appears as a dynamically 
generated resonance due to the strong interaction within 
and between the ttN and aN channels. This hadronic 
model has been quite successful in reproducing the t:N 
partial-wave amplitudes with total angular momentum 
J = 1/2 and 3/2 up to a cm. energy of 1.9 GeV. Recently, 
also the pole positions and residues of the resonances in 
this model have been extracted [54,, |55| . 

As alluded to above, in Ref. [i^ a dynamical coupled- 
channels model for pseudoscalar meson photoproduction 
based on the field-theory approach of Haberzettl [56J 
was introduced in conjunction with the Jiilich hadronic 
coupled-channels model . For the present application, 
however, we will use the reformulation of Haberzettl, 
Huang, and Nakayama i.50'] , which differs from the origi- 
nal approach of Haberzettl, Nakayama, and Krewald [4^ 
in some essential aspects that provide several practical 
advantages, as we shall explain in more detail in Sec. [Hi 

In general, the present photoproduction approach is 
distinguished from most existing dynamical models by 
the fact that it satisfies the generalized Ward-Takahashi 
identity for the production current [H, |57|] that ensures 
its full gauge invariance. By contrast, the vast majority 
of existing dynamical models at best provide only a con- 
served current but are not truly gauge invariant in their 
internal dynamics (see subsequent paragraph) . An ex- 
ception to this is the model of Refs. fSS-w], where gauge 
invariance is implemented following the prescription of 
Gross and Riska [6l|. See also Ref. [62|, where the issue 
of gauge invariance in pion photoproduction is discussed. 
An alternative method for achieving gauge invariance was 
developed in Ref. [G^I; however, it has not been applied 
in practical calculations so far. 



The present approach [i^, [sO, HI] is based on a fully 
microscopic (i.e., local) implementation of gauge invari- 
ance in which each electromagnetic current contribution 
associated with an internal subprocess of the reaction sat- 
isfies its own individual off-shell Ward-Takahashi iden- 
tity (WTI) — either as an ordinary WTI for a single- 
particle current |64| or as a generalized WTI for an in- 
teraction current [56l [STj — thus ensuring the overall 
gauge invariance of the physical current matrix elements. 
It was emphasized in Ref. [s^ that this feature is es- 
sential for making each individual contribution a con- 
sistent building block for the correct description of the 
reaction dynamics of the entire process. Most impor- 
tantly, these n-point-current building blocks (and their 
associated WTIs) remain the same whatever the dynam- 
ical context in which they appear. This is not simply a 
purely theoretical issue concerning the aesthetics of the 
formulation, but this consistency requirement across all 
possible reactions has immediate practical consequences. 
This was most clearly demonstrated in the recent study 
of the NN bremsstrahlung reaction [gl, |6^ whose fully 
gauge-invariant reaction amplitude was obtained within 
the same field-theory framework [5^ that provides the 
basis for the present formulation of pion photoproduc- 
tion ■ Specifically, it was shown in Refs. [GSl l66j 
that essential aspects of the bremsstrahlung process can 
be understood as time-reversed meson photoproduction 
processes NM -> N'^ (where M stands for a meson) 
whose dynamical details, therefore, can be described 
by the current building blocks from the correspond- 
ing meson-production formulation. The most important 
building block in this context turned out to be the four- 
point interaction current for the subprocess NM N'^ 
whose gauge-invariant construction in terms of its corre- 
sponding interaction-current WTI provided the contribu- 
tion necessary to resolve the longstanding discrepanc y of 
nearly a decade between the high-precision KVI data \&Q 
and the then existing models of the NN bremsstrahlung 
reaction. 

In the present work, going beyond the feasibility study 
given in Ref. (49| . we carry out an extended and more 
quantitative calculation of the jp — !■ 7r+n, 7p — !■ tt'^p, 
771 — ^ TT~p and jn — > Tr^n reactions using the input of 
the Jiilich ttN model of Ref. [33J to describe the hadronic 
final-state interactions (FSI) of the reactions. The pion 
photoproduction current itself is calculated using the 
novel formulation of Ref. (50i] whose details are given in 
the subsequent section. To keep the equations manage- 
able numerically, as a first step toward a more complete 
calculation, we explicitly include only the rjN and ttA 
channels and the basic ttN channel in the FSI loop inte- 
gration. Of course, the corresponding meson-baryon to 
meson-baryon T-matrices contain the information of all 
hadronic channels that comprise the Jiilich model. We 
calculate cross sections as well as beam asymmetries up 
to the total cm. energy ol W — 1.65 GeV. 

The present paper is organized as follows. In Sec.|lTl we 
briefly present our covariant formalism for pion photopro- 
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duction. In Sec. IIIIl we explain how the Jiilich hadronic 
model — which is based on TOPT — has been reformu- 
lated and matched to the present covariant formalism. 
In Sec. IIVI details of additional approximations made in 
our photoproduction model are presented. Our results 
for both the neutral and charged pion photoproduction 
reactions are presented in Sec. |Vl where some discussion 
is presented as well. Section IVII contributes to the dis- 
cussion about the theoretical uncertainties of our results. 
Finally, the summary and discussion of future develop- 
ments are given in Sec. I VIII Some details of the present 
model (interaction Lagrangians, form factors and the cor- 
responding parameter values) are given in the Appendix. 



II. FORMALISM 

The basic structure of pion photoproduction current 
M^^ seems fairly simple topologically. As shown in Fig.[l] 
it is comprised of the three s-, u-, and i-channel contribu- 
tions M^, M^, and , respectively, where the photon 
is attached to the external legs of the basic ttNN vertex 
and one diagram where the photon interacts with the 
interior of the vertex (correspondingly called the interac- 
tion current M^^^), i.e.. 



(1) 



In general, all four current contributions here are fully 
dressed. Following Ref. [gl], the first three (pole-type) 
diagrams are usually called class-A diagrams and the last 
(non-pole) diagram is of class B. This simple topological 
structure is also reflected by the corresponding basic tree- 
level Feynman diagrams for the process when reducing 
the full complexity of the current to using bare prop- 
agators and vertices only. In full detail, however, the 
microscopic dynamics of the production process is much, 
much more complex because of dressing effects of the 
hadron propagators and vertices. Moreover, since parti- 
cle number is not conserved (e.g., internally any number 
of mesons can partake in the process) , the full process is 
highly non-linear as a matter of course. 

Comprehensive theoretical formulations of photopro- 
duction processes must be able to incorporate these dy- 
namical complications at least in principle. In addition, 
the corresponding production current Af^ must obey 
gauge invariance as the manifestation of U(l) symmetry 
which is of fundamental importance for any photopro- 
cess because it provides a conserved current and thus 
implies charge conservation. For the pion production 
current , in particular, gauge invariance is formu- 
lated in terms of the generalized Ward-Takahashi identity 




Figure 1. Generic structure of the pion photoproduction cur- 
rent for 7A'' — >■ TiN according to Eq. ([T]). Time proceeds 
from right to left. Nucleons and pions are depicted by solid 
and dashed lines, respectively, and the photon is shown as 
a wavy line. The four-momenta shown at the external lines 
are those of the respective particles used in this work. The 
first three diagrams comprising class-A diagrams ^60| depict 
the S-, U-, and t-channel pole diagrams , M^, and , 
respectively, where s, u, and t are the respective Mandelstam 
variables of the internally exchanged particles, as indicated. 
The last diagram (called class B [68]) shows the contact-type 
interaction current Mlt^. 



where the four-momenta are those shown in Fig. [TJ The 
vertices here correspond to the (fully dressed) ttNN 
vertices F in the respective kinematic situations corre- 
sponding to the Mandelstam variables x = s,u, t, as 
shown in Fig. [TJ The propagators for the nucleon and 
pion are denoted by S and A^, respectively, and the 
charge operators for the initial and final nucleon and for 
the outgoing pion are Qi, Q/, and Q^, respectively. Ob- 
viously, this expression vanishes for on-shell hadrons and 
thus provides a conserved current. This off-shell formu- 
lation of gauge invariance, however, goes beyond that by 
providing a local constraint on the gauge invariance of 
the photoproduction current that is similar to requiring 
the usual Ward-Takahashi identity for the single-particle 
currents 6j|, which for the nucleon reads 



k^J^' = S-\p + k)QN 



QnS-\p) , 



(3) 



where is the nucleon's generic charge operator. Both 
requirements © and ^ (and its analog for the pion) 
are essential for the internal consistency of microscopic 
formulations of photoprocessesQ 

One can easily show that assuming the validity of the 
usual WTI ([3]) for the electromagnetic nucleon current 
Jf^ and its analog for the pion current, the generalized 
WTI ^ implies 



^Ct = -P^e^ + Fuef + Fte^ , 



(4) 



where the F^ are the vertices F^ of ([2]) stripped of their 
isospin operators r that now appear in 



Cj = tC 



ef = QfT, and e^r = QttT (5) 



which are the charges for all external hadron legs in an 
appropriate isospin basis (with all corresponding indices 



-F,S{p + k)Q,S~\p) 
+ S-\p')QfS{p' ~k)Fu 
+ A-\q)Q^A^iq-k)Ft 



(2) 



^ As mentioned in the Introduction, a particularly striking exam- 
ple of this, with immediate practical consequences for the de- 
scription of experimental data, was recently discussed for NN 
bremsstrahlung | 65t . ,66] . 



and summations suppressed). In other words, the rela- 
tion 

ei=ef + (6) 

describes charge conservation for the pion photoproduc- 
tion process. With the single-particle WTIs for nucleons 
and pions given, Eq. ^ is completely equivalent to the 
generalized WTI ©. It is this Eq. g]) for the contact- 
type interaction current M/^^, in particular, that is being 
exploited here for the purposes of preserving the overall 
gauge invariance of the production current. 



A. Dyson-Schwinger framework 

The complete (i.e., fully dressed, non-linear, and 
gauge-invariant) structure of pion photoproduction was 
described by Haberzettl [s^ within a covariant field- 
theoretical Dyson-Schwinger framework. In practice, 
however, the complexity of the full formalism needs to be 
truncated at some level to make it numerically tractable. 
Since doing so invariably leads to a violation of gauge in- 
variance, one must find prescriptions to restore this fun- 
damental symmetry. And one should do so in a manner 
that incorporates as many of the original reaction mech- 
anisms as possible. Particularly relevant in this context 
are the hadronic final-state interactions of the outgoing 
meson-nucleon states since these form essential parts of 
the dynamical content of the interaction current Af/^j 

[Hi. 

In Ref . [401 , it was shown how to approximate the full 
formalism in a manner that includes the full hadronic 
final-state interactions while at the same time preserv- 
ing the gauge invariance by reproducing the generalized 
Ward-Takahashi identity (HJ for the production current. 
In the present paper, we will follow the variant of the pro- 
cedure of Ref. 149^] put forward recently in Ref. [13] . We 
emphasize that without any approximation, i.e., as full 
Dyson-Schwinger formulations, the respective results of 
Refs. [i^jls^lsg are completely equivalent. However, the 
features of Ref. [50] turn out to be particularly well suited 
for approximating the dressing effects of the electromag- 
netic nucleon current J'^ in a way that is reciprocally 
consistent with the production current itself while at 
the same time preserving local gauge invariance. Here, 
we will recapitulate the details of the approach given in 
Ref. [lOl only insofar as it is necessary for the present 
application. To this end, we make extensive use of dia- 
grammatic representations. Full details and derivations 
can be found in Ref. [13], and in Refs. [igUsj. 

Note that on the hadronic side, for simplicity, we will 
speak here most of the time explicitly only of pions and 
nucleons. However, these particles are to be taken as 
representatives for all mesons and baryons, respectively, 
that take part in a particular application. Correspond- 
ingly, any explicit equation appearing here for the t:N 
channel is to be read as a matrix equation that couples 
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Figure 2. Microscopic structure of the photoproduction pro- 
cess according to Ref. [sO]. The line styles are the same as 
in Fig. [T] with thick solid nucleon lines depicting the propa- 
gation of dressed intermediate nucleons. Solid and open cir- 
cles indicate dressed and bare vertices, respectively, (a) The 
full current with its FSI contribution written as a loop 
integration over the full nN T-matrix. (b) Details of the 
Born-type current of Eq. ((8]) that comprises the u and 
t-channel contributions and Mj^ , and the contact-type 
four-point current depicted as a solid square whose de- 
tails in turn are shown in part (c). The first term here is 
the KroU-Ruderman contact term rnKp. The various hadronic 
boxes labeled U depict the non-polar ttN irreducible driving 
term of the ttN T-matrix, as shown in Fig. [^c). The partic- 
ular five-point current C/^ appearing in the first loop arises 
from attaching a photon to U ; its lowest order is shown in 
part (d). The four-point currents appearing here correspond 
to the full interaction current M^^^^ of Fig. [T] The box labeled 
Bt in (a) and the loop integrations labeled L in (c) indicate re- 
strictions of the corresponding B'^ contributions to transverse 
and longitudinal pieces, respectively (see text). The details 
of the full nucleon current J** (solid circle) and its auxiliary 
s-channel restriction (solid square) are depicted in Fig. [S] 

all hadronic meson-baryon channels once the full com- 
plexity of the reaction dynamics is turned on. 

Following Ref. [13] , the full pion photoproduction cur- 
rent AI^ in the one-photon approximatior(3 can be writ- 
ten as 

M'" = FsSJ^ + + TGoBf^ , (7) 

with the non-pole Born-type current given by 

Bf" = + -I- . (8) 

These equations are depicted diagrammatically in Fig. [21 
The (fully dressed) u- and i-channel contributions AII^ 



^ This is not a serious limitation since higher-order contributions 
are suppressed by ~ 1/137. For most if not all practical 
purposes, therefore, the one-photon approximation is perfectly 
adequate. 
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(b) 



Figure 3. (a) The dressed electromagnetic nucleon current J*^, with the first term on the right-hand side (depicted as a solid- 
square three-point current) corresponding to appearing in Eq. [T] (b) The details of Jl^ . The first term (with open circle) 
is the bare current and the second term contains the KroU-Ruderman term (open-circle four-point current) integrated over by 
the dressed nNN vertex (solid circle; cf. Fig. [5]). The remaining loop contributions only contribute longitudinally (as indicated 
by the index L). By contrast, the loops in (a) only contribute transversally, as indicated by the letter T. The solid-square 
contact-type four-point current in the loops at the end of both lines is given in Fig. [2jc). 



and Mj^ of Eq. ([T]) appear here unchanged, but the s- 
channel current and the interaction current Afj^^ are 
now spread out over the remaining terms. The first term, 
FsSJ!^, contains part of the s-channel pole contribution 
via the fully dressed nucleon propagator S, however, with 
an electromagnetic current contribution for the nu- 
cleon that is only part of the full nucleon current J^. 
The situation is depicted in Fig.[3]and will be discussed in 
more detail below. The contact term [see Fig. ^c)] 
contains part of the full interaction current M^^^. The 
remaining pieces for both and M-^^. come from the 
TGqB!^ term which describes the ttN final-state interac- 
tion mediated by the loop integration over the full t:N 
T-matrix, where Go is the intermediate propagator of the 
free ttN pair within the loop. Note that only the trans- 
verse parts of the contributions of these loop integrations 
are to be taken into account here, as indicated by the in- 
dex T0 The respective missing pieces for and M^^^ 
arise from splitting the full T into an s-channel pole part 
and a non-pole part according to 

T=\F)S{F\+X, (9) 

where the first term contains the s-channel nucleon pole 
via the dressed nucleon propagator S and X denotes the 
remaining non-pole part|f| This splitting is depicted dia- 
grammatically in Fig. 2] together with the Bethe-Salpeter 



^ For definiteness, since transverse parts are not unique, through- 
out this work here the index T on a current indicates j!^ = 
s^Evj" je^ where £^ is the (transverse) photon polarization vec- 
tor. The splitting into longitudinal and transverse parts satisfies 

* We follow here the notation of Ref. [5(1 . In other words, for 
notational simplicity, we do not employ the often-used notation 
and for the pole and non-pole contributions of T, re- 

spectively, since this notation tends to make the equations dif- 
ficult to parse. For similar reasons, to avoid the extensive use 
of adjoint daggers, we use the ket and bra notation \F) and (F 
to denote (dressed) nNN vertices that describe — ^ ttN and 
ttN — > N, respectively. Within the context of Eq. (0, therefore, 
the s-channel vertex Fs may also be written as |F). The bras 
and kets, however, are not to be misconstrued as Hilbert-space 
vectors. 




Figure 4. (a) Representing the full T-matrix for nN scattering 
via its (fully dressed) s-channel pole part and its non-polar 
contribution X, according to Eq. ((9]). (b) The Bethe-Salpeter 
scattering equation for the non-pole amplitude X, Eq. (|10[) . 
Its driving term U comprises all ttN irreducible mechanisms. 
The corresponding lowest-order contribution is the M-channel 
exchange shown explicitly in (c). Attaching a photon to this 
diagram provides the five-point current given in Fig. [SJ^d). 

equation for the non-pole amplitude, 

X = U + UGqX, (10) 

driven by the (fully dressed) non-pole nN irreducible 
mechanisms subsumed in U whose lowest-order contribu- 
tion is the It-channel exchange shown in Fig. Iljc) . The 
_F's here are the fully dressed irNN vertices related to 
the bare vertex / by 

\F) = \f)+XGo\f) . (11) 

This dressed vertex is represented diagrammatically in 
Fig. [5] together with the dressing mechanism for the nu- 
cleon propagator 

S^So + SoJ:S (12) 

given in terms of the bare propagator 5*0 and the nucleon 
self-energy loop S. 

The entire set of equations summarized by Figs. [2] 
through [5] provides a complete solution of the problem by 
forming a tower of coupled non-linear Dyson-Schwinger- 
type equations that can only be solved iteratively. View- 
ing this coupled set in its entirety for a given bare input. 
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Figure 5. Dressing mechanisms for (a) the nucleon propagator 
S in terms of the bare propagator So (thin line) and the self- 
energy loop and (b) the ttNN vertex F according to Eq. . 
Dressed and undressed vertices are shown as solid and open 
circles, respectively. 

the splitting © of T into pole and non-pole contributions 
is unique. However, this uniqueness generally is lost once 
approximations are introduced. As a consequence, the re- 
sulting individual pieces of a particular approximation of 
the splitting ^ may exhibit undesirable numerical arti- 
facts which the T matrix itself does not possess. (This is 
the case, for example, for pole and non-pole parts result- 
ing from the Jiilich model we use here [53] •) In order to 
avoid such problems to some extent one may, instead 
of Eq. ([T0|) . employ the Bethe-Salpeter equation 

T = V + VGoT (13) 

for the entire T-matrix, with the driving term 

V^\f)So{f\ + U , (14) 

where / and Sq are the bare vertex and bare propagator 
of Eqs. ([TT|) and respectively. 

The photoproduction current of Eq. ([7]) is formulated 
assuming that the hadronic scattering information en- 
ters the problem in terms of T, i.e., in terms of a ma- 
trix coupled-channels equation of the generic type given 
in (1131). An equivalent formulation in terms of X e:^- 
ists 'so*, 'sg'I , however, since the usual hadronic coupled- 
channels approaches (like the Jiilich model we employ 
here) produce full T-matrices as a matter of course, it 
is of considerable practical advantage to utilize this in- 
formation directly for the photoproduction current. (See 
also the corresponding discussion at the end of Sec. IIIII ) 



We mention that in the one-photon approximation, one can- 
not formulate the photoproduction process without employing 
a sphtting of T into pole and non-pole contributions of some 
sort. The two opposite extremes of such splittings are, at one 
end, the Dyson-Schwinger-type splitting ^ and, at the other 
end, the splitting where the pole part only contains the phys- 
ical pole position itself, without any dressing effects, and the 
(constant) residue at this pole, with everything else appearing in 
the non-pole part. Going beyond the one-photon approximation, 
one may perhaps be able to avoid the aforementioned problems 
if one makes the -yN channel an integral part of the set of cou- 
pled channels thus producing the photoproduction current as a 
solution element of T when solving the corresponding coupled- 
channels matrix equation jTsjl. However, we have no proof of 
this conjecture. Moreover, preserving gauge invariance in that 
approach is much more complicated than in the present one. 



B. Gauge-invariant truncation of the pion 
photoproduction current 

As mentioned, the preceding formulation is exact, but 
it is also very complex in detail and therefore, at present 
at least, cannot be solved in its entirety without making 
some approximations to render the problem manageable 
numerically. In the particular formulation {SQ] used here, 
there are two obvious places for such approximations — 
the dressed nucleon current depicted in Fig. [3] and the 
contact-type five-point current Af^ shown in Fig.[2][c). 

The detailed (exact) structure of the electromagnetic 
nucleon is given by [13] 

= Jl^ + {F\ Go {M^ + Mt + Ml^X , (15) 

as one easily reads off Fig. [3] The currents within the 
loop structure here are given entirely by pieces occurring 
in the pion production current itself. As was shown in 
Ref. 50(] , this reciprocal consistency between the nucleon 
current and the production current was essential 
for deriving the particular form of given in Eq. ([7]) 
because the transverse loop contributions of Eq. (|15l) . 
in particular, survive unchanged in the loop integration 
that provides the FSI contribution in Eq. ([7]). This form, 
moreover, easily permits to treat both and with 
a consistent set of approximations. 

1. Approximating the nucleon current 

First, we note that all loop integrations in Eqs. ([7]) 
and (|15p contain the full nucleon current within the u- 
channel loops. The nucleon current, thus, couples back 
into itself and these contributions, therefore, are part of 
the non-linear hierarchy of Dyson- Schwinger equations 
that, in principle, needs to be solved iteratively requir- 
ing prohibitive amounts of computational resources. To 
avoid this complication, we truncate these loop contri- 
butions by employing in our present application only the 
usual (on-shell) and a^'^ki, pieces for the nucleon cur- 
rent within all the loops. This approximation step thus 
linearizes the nucleon-current contributions. Second, we 
note that the current piece that appears in the s- 
channel part of ([7]) obeys the same WTI ([3]) as the full 
current, i.e., 

fc^J^ - kf^J^ - S-\p + k)QN - QnS-\p) , (16) 

since the their difference in ([T5|) is purely transverse. The 
detailed expressions for J^, moreover, mainly are com- 
prised of longitudinal contributions that arise from the 
photon being attached to the self-energy loop of the nu- 
cleon propagator [cf. Figs.[3Jb) andlSfa)]. 

To reproduce the WTI (fTB)) . it is convenient to utilize 
the Ball-Chiu current [69j . 

J, ^ ,^,u S-'{P + k)QN-QNS-Hp) 
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(2p + fc)^ 

s — 



Qn 



A{s) + A{p^) 



(17) 



where s = (p + fc)^ and A is one of two scalar dressing 
functions {B being the other) resulting from the generic 
form 



S{p) 



1 



;5A(p2) - jnB{p'^) 



(18) 



of the dressed nuclcon propagator; m is the nucleon mass. 
The dressing functions A^p^) and B{p'^) are constrained 
to produce a unit residue at the nucleon pole where p^ — 
[H, nil . By construction, is the minimal current 
that satisfies the WTI ((T6)) for fully dressed propagators 
and it is both non-singular and symmetric. Without lack 
of generality, therefore, we may write 



J<: = Jf + Tt: 



(19) 



where is the transverse remainder defined by this re- 
lation. 

Within the s-channel context of Eq. ([7]), the Ball-Chiu 
current provides a particularly simple expression when 
considering the half-on-shell situation, with an incoming 
nucleon spinor u{p) on the right and an outgoing propa- 
gator S{p -\- k) on the left. One easily finds 



SJ!fu = 



1 



y + 



K2 - Ki 



2m 



]nu{p), (20) 



where p^ 



Albeit written in a somewhat different 



way, this expression is equivalent to Eq. (17) of Ref. [50|, 
and the two independent coefficient functions Hi = Ki{s) 
{i = 1,2) that are directly related to the propagator 
dressing functions A and B are given in Eq. (19) of 
Ref. [50] ■ We omit their details because we will not make 
use of them here. We only mention that the on-shell 
values at s = w? for both coefficients are identical, i.e., 
ki(to^) = K2{m^) = A(m^) — 1, and that one easily finds 
that {ii2 — Ki)/{s — Tin?) possesses a finite limit for s ^ 
if we assume that the dressing functions A and B are ana- 



lytic functions in the vicinity of s 



This means they 



both vanish in the structureless limit [where A{m?) = 1], 
thus leaving in (|20p only the usual 7'^ Dirac current to- 
gether with a structureless propagator. All efi^ects of the 
dressing thus reside in the terms that depend on the Kj 
whose overall contributions are manifestly transverse. In 
the present application, we will absorb the (z = 1,2) 
dependence in some fit parameters. 

We note in this context that the Ball-Chiu current J^" 
does not fully contain the anomalous- moment contribu- 
tion of the usual Pauli part of the on-shell nucleon cur- 
rent. These anomalous contributions arise from other 
transverse dressing mechanisms. Writing the on-shell 



matrix element of the full nucleon current between nu- 
cleon spinors in the usual manner as 



N 



2m 



(21) 



where e is the fundamental charge unit, kjv is the anoma- 
lous moment of the nucleon, and (5 at — 1,0 for the pro- 
ton and neutron, respectively, the anomalous part can be 
written as 

eu I — ; kn \ u — u \ eOM — '«o 



2m 



2m 



+ {F\ GoB^ 



(22) 



where kq = Ki(m^) = K2(™^) — A(to^) — 1. This result 
follows from Eq. (|15l) utilizing Eq. ([T5|) and the expression 
performing expansions of Ki{s) around s = m^. In 
principle, this equation is exact if T'^ and B!^ could be 
calculated without approximations. We see here that for 
the proton, the Ball-Chiu term ko contributes partially 
to the anomalous moment, but for the neutron, it does 
not contribute at all. 

In Ref. [131 , it was advocated to put = and ensure 
the anomalous contributions by adding a contact term to 
the approximation of to be discussed in the follow- 
ing subsection in the context of Eq. (l30l) . In the present 
application, however, we proceed differently. We expand 
the full dressing functions Kt appearing in the Ball-Chiu 
contribution (j20p around their on-shell points and use the 
resulting coefficients as fit parameters. In addition, we 
assume that we can approximate by a single trans- 
verse contact term whose operator structure is given by 
a^^k^ alone. It is the corresponding coefficient, in partic- 
ular, that ensures that we can reproduce the anomalous 
moments. For the half-on-shell matrix element of ap- 
pearing in the s-channel pole term of Eq. we then 
obtain the approximation 



SJ!!u SJi'u = e 



^a^"'k^ 



2m 



tt'^k^ 
4m^ 



.a'^'k^ 
4m^ 



C2 



u{p) 



where 



7'^(s-m2) - {2p + kYlf: 
2m 



(23) 



(24) 



is a transverse current operator that follows from the 
expansions of the Ball-Chiu contributions. For each nu- 
cleon, this approximation contains three dimensionless 
coefficients: a real parameter kq, and two complex num- 
bers Ci and C2 for the contact contributions. 

We emphasize that the approximation just discussed 
only pertains to the true nucleon current. In the full 
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coupled-channels treatment, the structure of Eq. (O also 
comprises contributions from currents Jjf that describe 
nucleon-resonance transitions. Such transition currents, 
however, must necessarily be transverse and thus corre- 
spond only to the first two diagrams on the right-hand 
side of Fig. [3][b). In the present applications, thus, we 
describe these transition currents by the point vertices 
of the nucleon-resonance photo-transition Lagrangians 
given in Eqs. (jA.lSp and (jA.16|) of the Appendix. 

We also emphasize that the present approximate treat- 
ment of the s-channel term (|23p that permits one to 
describe the corresponding fully dressed contribution in 
terms of three parameters requires the complete recip- 
rocal consistency between the nucleon current and the 
photoproduction current, as derived in Ref. [s^l- The 
approximation scheme of Ref. 49], by contrast, is am- 
biguous in its treatment of certain transverse current 
contributions that have no bearing on gauge invariance. 
Specifically, the undetermined transverse current ap- 
pearing in Eq. (21) of Ref. [1^ was taken to be zero in 
the preliminary applications reported there. To achieve 
complete structural equivalence with the present results, 
one finds that T'" must be chosen as T^' = UGo{M^)^ in- 
stead if one wants to preserve the equivalence even when 
making approximations. 

2. Approximating the four-point contact current 

In view of the fact that the FSI loop in ([7]) is trans- 
verse and that the approximation of the s-channel current 
J^" by the Ball-Chiu current leaves the correspond- 
ing WTI unchanged, the entire burden for reproducing 
the generalized WTI for the production current rests 
with the contact-type four point current whose four- 
divergence, by construction, must be the same as that of 
the interaction current M^^^ in Eq. i.e., 

K ^Ct = K^^c = ~Fse^ +Fuef+Fte^. (25) 

As far as preserving gauge invariance is concerned, there- 
fore, the situation is very much the same as in the ap- 
proach of Ref. , even though the details of the present 
formulation are somewhat different. In other words, any 
approximation of must satisfy the constraint (l?5|) . 

The internal dynamical mechanisms of as depicted 
in Fig. [2jc) are quite involved and cannot, at present, be 
incorporated fully in numerical applications. In Ref. [49j, 
it was discussed in detail how to find approximations of 
the mechanisms within at any desired level of so- 
phistication while still reproducing ((25|) . For the present 
application, we simply quote the results of [49{ valid 
when the t:NN vertices are described by phenomeno- 
logical form factors, as it is the case here. 

Phenomenologically, the ttNN vertex stripped of its 
isospin dependence can be written as 

Fx = gNNTv-!r^fx , (26) 
2m 



where we take here pure pseudovector coupling, with 
gNNrr being the physical coupling constant. The sub- 
script X — s, u, or t indicates the kinematic context; 
is the phenomenological form factor for the correspond- 
ing reaction channel normalized to unity when all hadron 
legs are on-shell and q is the four-momentum of the out- 
going meson. Following Ref. 143], we define an auxiliary 
current 

t — u — p''^ 

- e.^^^{2p + kr , (27) 
s — p'' 

with 

F^l~hil- SJs) (1 - Sufu) (1 - Stft) , (28) 

where the four-momenta are those given in Fig. [TJ The 
constant is unity if the corresponding x-channel con- 
tributes to the reaction in question, and zero otherwise. 
The parameter h may be an arbitrary (complex) function, 
h = h(s, u, t), possibly subject to crossing-symmetry con- 
straints. In the present work, it is simply taken as a fit 
constant for the sake of simplicity. The four-divergence 
of the auxiliary current evaluates to 

k^Cf^^^fse. + Uef + fte^. (29) 

We emphasize that is non-singular since, by construc- 
tion, the propagator singularities are canceled by the cor- 
responding zeros of fx — F. In other words, C'^ is a true 
contact current. 

The gauge-invariance preserving (GIF) approximation 
of the contact current Af^ then may be chosen as [igllsoj 

M^^^MH^g^^J-^^i^C^ 
2m 

-gNNn^{enft-f3kpCP) , (30) 
2m 

where /? is a free fit parameter. One can easily check that 
this choice for satisfies the gauge-invariance condi- 
tion Essential in this respect is the relation (|29l) . 
Note that the term proportional to the pion charge e^r in 
(pO)) is a dressed version of the familiar Kroll-Ruderman 
contact current where the dressing form factor ft is that 
of the i-channel amplitude Mj^. The original (undressed) 
Kroll-Ruderman term survives in this GIF current since 
ft — ^ + {ft — 1)7 i-e., the dressing is given by the addi- 
tional ft — ^ contribution. 

We emphasize that the present formulation is equally 
valid for real and virtual photons. For pion photopro- 
duction with real photons, in particular, only transverse 
currents contribute to the physical amplitude. The lon- 
gitudinal loop contributions shown for Af^ in Fig. (^fc) , 
therefore, are projected out and only the (bare) Kroll- 
Ruderman current and the loop over the the five-point 
interaction current [/^ of Fig. [IJd) remain. Hence, sub- 
tracting the usual (bare) Kroll-Ruderman term from the 
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GIP current (|30p . the remaining expression is a phe- 
nomenological approximation of the latter loop contri- 
bution for real photons. 

Finally, we mention that when we consider interme- 
diate channels other than ttN (cf. Sec. IIVI below), the 
corresponding contact currents Mj^ that appear in the 
transverse loops TGqB!^ of Eq. ([7]) are easily constructed 
in analogy to the explicit TriV case given here by using 
the generic constraint (p5)) . 



3. Covariant three-dimensional reduction 

The present formalism is a fully covariant approach. 
To calculate any reaction amplitude in a full four- 
dimensional framework is a daunting numerical task. 
Therefore, in the present work, we approximate the full 
meson-baryon two-body propagator Go by the Kady- 
shevsky propagator ^7Q|], as specified in Eqs. and ([57)1 
in the next section. This propagator, together with an 
implied energy delta function, restricts the propagation 
of the intermediate meson and baryon to their respective 
energy shells, thus reducing the four-dimensional loop in- 
tegration to a three-dimensional one without destroying 
the covariance of the equation. This reduction is also nec- 
essary to be able to match up the hadronic Jiilich-model 
input with the present formalism. 



III. MAKING THE JULICH nN MODEL 
COMPATIBLE WITH THE PRESENT 
APPROACH 

In the present work, we employ the Jiilich dynami- 
cal coupled-channels nN model |3ll - [33| to account for 
the hadron dynamics. As mentioned in the Introduc- 
tion, this model is formulated within time-ordered per- 
turbation theory (TOPT) which is a non-covariant three- 
dimensional formalism [5l|. We, therefore, need to pro- 
vide a procedure to match it to the covariance of the 
present photoproduction formalism. 

To this end, we first note that the Jiilich T-matrix 
satisfies the three-dimensional Lippmann-Schwinger-type 
equation 

7to(p',P;^) = Vto{p',p;z) 

+ J dV'^To(p',p";^)GTo(p",^)TTo(p",p;^) , 

(31) 

where z denotes the total-energy variable and TO indi- 
cates that all entities result from the time-ordered formal- 
ism. The intermediate pion-nucleon propagator reads 



Gto{p",z) 



1 



(32) 



z- E{p")-uj{p") + iO 
- p"2 stands for the energy of the 



where E{p") = ^ m?- 
nucleon and i^(p") = \J m\ + p"^ for the energy of the 



Table I. Cutoff parameters (in MeV) in the form factors of 
Eqs. (fOTll and (|XT8)) . 



A 



A, 



A„ 



600 



725 



1343 



1645 



Table II. Effective electromagnetic coupling constants for 
isospin 1/2 resonances. The subscripts p and n denote the 
resonance coupling to 7p and 771, respectively. 





Di3(1520) 


5'ii(1535) 


Sii(1650) 


Pi3(1720) 




-1.23 


-1.14 


0.02 


4.23 




1.58 






-12.56 


yn 


2.69 


0.21 


-0.17 


-0.73 


(2) 
9n 


-3.58 






2.52 



Table III. Effective electromagnetic coupling constants for 
isospin 3/2 resonances. 

P33(1232) 53i(1620) D33(1700) P3i(1910) 



2.01 

4.48 



1.63 



-6.56 
5.38 



-1.44 



pion, with the mass m^. The complex notation iO indi- 
cates the physical limit to the upper edge of the scatter- 
ing cut to provide proper boundary conditions. Equa- 
tion pi|) written in the cm. system is structurally very 
similar to the Bethe-Salpeter equation (fT3| after it has 
been subjected to a covariant three-dimensional reduc- 
tion. Therefore, to match the two, one needs to ensure 
that Eq. pip transforms covariantly away from the cm. 
system. 

We define now 



V{p',p;V^s) 
Tip',p;^s) 



(27r)V(p')^To(p',P;^)p(p) , (33) 
{2nfp{p')TToip\p;z)pip), (34) 



where 



p(p) = 2/b(pHp) 



(35) 



corrects the non-covariant normalizations of the TOPT 
plane-wave states and where the invariant mass ^/s is 
defined as ^/s = z in the cm. system. Equation pip 
then may be recast in the equivalent form 

r(p',p;V^) =F(p',p;V^) 



(27r)3 



V{p', p"; V^) Go(p", V~s) T(p", p; ^) 



where 



Go(p",V^) = 



p2(p") ^- E{p")-uj{p")+iO 



(36) 



(37) 
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Table IV. Dimensionless parameters Ci, C2, and kq in the 
electromagnetic nucleon current, Eq. (I23|l . 





C2 




proton 2.07 e"^-^*^ 
neutron 12.35 e^'^^* 


34.62 e° °®' 
25.96 e^ ™' 


-4.24 
3.24 


Table V. Dimensionless parameters h and /3 in the contact 
current , Eq. (|30p. for the various meson-baryon channels. 




TT J3 TIP 


ttA 


/i 1.43 1.44 
P -0.17 -1.43 


2.01 -10.04 
-0.10 3.45 


0.74 
-1.75 



Equation (136]) is formally identical to the Kadyshevsky 
equation [70J which is the result of a particular covari- 
ant three-dimensional reduction of the Bethe-Salpeter 
equation that preserves elastic unitarity. Hence, utilizing 
the time-ordered Jiilich T-matrix in the redefined form 
P4p is completely consistent with the covariant three- 
dimensional reduction of the photoproduction current 
given in Eq. ([7]). For the corresponding reduced meson- 
baryon propagator one must use the form (j37p . of course. 

We recall in this context that pions and nucleons are 
used here as simplified generic tags for all mesons and 
baryons incorporated in the Jiilich model. Therefore, 
all entities discussed here need to be considered as ele- 
ments of matrices labeled by meson-baryon channels and 
all equations contain implied summations over all ele- 
ments of this channel space. Moreover, for quasi-two- 
body channels criV, pN , and ttA, the propagator ([57]) is 
to be slightly modified by including the self-energy of the 
corresponding quasi- free particles [31, 32] • 

It was alluded to above in the context of the Bethe- 
Salpeter equation ([T3l) that the splitting of T into pole 
and non-pole contributions is no longer unique if the re- 
spective pieces are evaluated in a framework that trun- 
cates the non-linearities. Owing to this arbitrariness, the 
pole and non-pole parts of the Jiilich model both exhibit 
spikes in the Pn partial wave near the ttN threshold due 
to the presence of a nearby unphysical pole in the non- 
pole transition matrix X. This pole in X also affects 
the pole part of T since the dressed vertices F appear- 
ing in the splitting (|9]) also depend on X via Eq. ([TT|). 
and one finds that in the sum of pole and non-pole parts, 
these spikes cancel each other precisely to yield a smooth 
full T-matrix [54i] . The spikes, therefore, are unphysi- 
cal artifacts of the particular way the pole and non-pole 
contributions are calculated in the Jiilich model. To cir- 
cumvent problems of this kind, we have chosen here, in 
Eq. ([7]), to work with the full T-matrix to describe the 
hadronic final-state interaction. This, however, does not 
eliminate completely all numerical difficulties because the 
s-channel contribution in Eq. ([7]) still involves the dressed 
NNtt vertex which, as explained, also exhibits the 




(b) 



Figure 6. Hadronic states taken into account in the present 
model (see Sec. IIVI for details) : (a) intermediate s-channel 
baryons for 7A'^ — >■ ttN {R subsumes all resonances except the 
A) contributing to Fig. [2ja); u- and f-channel contributions 
as they appear in B'^ of Fig. [2{b) for (b) 7A'' — > tvN, (c) 
7A'^ — >■ r/N, and (d) 7A'^ — >■ ttA. Diagrams (c) and (d) only 
contribute within the FSI loop. 



spike by itself. To avoid this problem, we replace here the 
dressed vertex Fg of the Jiilich model by the NNtt vertex 
obtained within the Feynman prescription for the effec- 
tive NNtt interaction Lagrangian /^ArATTr given in the Ap- 
pendix employing the physical NNtt coupling constant 
and physical nucleon mass. 



IV. MESON-BARYON CHANNELS INCLUDED 
IN THE PRESENT MODEL 

In the present work, as a first step towards a more 
complete calculation, only the ttN, rjN and ttA channels 
are included as intermediate states in the loop integral 
in Eq. ([7]). However, this restriction only concerns the 
photoproduction sector. In the hadronic sector, all of 
the five channels ttN, rjN, ttA, aN, and pN are included 
via the channels that couple into all of the Jiilich T- 
matrices. The main reason we do not include the aN and 
pN channels of the hadronic Jiilich model in the present 
photoproduction calculation is purely practical, namely, 
the corresponding transition couplings of the photon to 
these channels would require additional free parameters 
beyond what we can handle numerically at the moment. 
Moreover, for the pN channel, we expect significant con- 
tributions only at energies higher than the cm. energy 
oi W = 1.65 GeV considered in the present work. 

The channels and intermediate hadrons that con- 
tribute explicitly in the present calculation are shown 
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Figure 7. (Color online) Differential cross section for 7p — ^ n^n as a function of the scattering angle. The first and second 
number in each pair of parentheses correspond to the photon laboratory momentum and the nN cm. energy, respectively. 
Data are taken from Refs. [73-[74ll. 



in Fig. |6l As depicted in Fig. EJa), in addition to 
the nucleon that is treated as described in Sec. IIIB 1[ 
there are eight resonant contributions that enter the 
s-channel term of Eq. ©, i.e., 5'ii(1535), S'ii(1650), 
53i(1620), P3i(1910), Pi3(1720), i?i3(1520), ^33(1232), 
and 1)33(1700). All of the corresponding hadronic ver- 
tices Fs and the dressed resonance propagators S are 
taken from the Jiilich TriV model. The correspond- 
ing electromagnetic nucleon-resonance transition vertices 
should be calculated from the analog of Fig. ^h). 
However, since the transition currents are all transverse 
by themselves, the analogs of the longitudinal loops in 
Fig.[3l^b) vanish. Moreover, since we do not consider bare 
KroU-Ruderman-type four-point currents for the reso- 
nances, we are only left with the bare currents that are 
given by the corresponding effective Lagrangians (jA.15[) 
and (jA.16p in the Appendix. This is one of the gratifying 
features of the present formulation that greatly simplifies 
the calculation of the s-channel resonant contributions. 

The ttN contributions of Fig. ^h) enter both 

the Born-type current and the FSI loop TGqB!^ of 
Eq. ([7]). The contributions of Fig. IHlJc) and Fig. [BUd) for 



7iV — > rjN and — > ttA, on the other hand, only en- 
ter the FSI loop. For the ^-channel currents, we include 
the contributions from the tt, p, uj, and ai exchanges in 
the — > ttN channel, from p and w exchanges in the 
7A'^ — >■ 7]N channel and from tt and p exchanges in the 
7iV — ttA channel. These diagrams are calculated using 
the corresponding Lagrangians listed in the Appendix. 
All the hadronic cou pling constants are taken from the 
Jiihch TTiV model [sJil [nl- The MM'7 electromag- 
netic coupling constants are determined from the radia- 
tive decay of the relevant mesons in conjunction with 
SU(3) symmetry considerations. The numerical values of 
the coupling constants are given the Appendix in connec- 
tion with their associated interaction Lagrangians. The 
only adjustable parameters here are the cutoff parame- 
ters in the off-shell form factors introduced at the elec- 
tromagnetic vertices. These form factors are also given 
in the Appendix. 

The only contribution to the u-channel amplitude A/^ 
considered in the present work is that involving the nu- 
cleon (nucleonic current) and A in the intermediate state. 
Contributions to the w-channel from the baryon reso- 
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Figure 8. (Color online) Differential cross section for 7p — >■ n^p as a function of the scattering angle. The first and second 
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Data are taken from Refs. [73471 . 



nances are not included for consistency reasons since the 
Jiilich hadronic model (and most other existing dynam- 
ical models for that matter) does not include the reso- 
nance contributions to the u-channel except the A(1232) 
resonance. Apart from an enormous numerical demand 
in keeping the self-consistency between the s- and u- 
channel amplitudes due to the dressed vertices and prop- 
agators, the u-channel contribution from the baryon reso- 
nances are not expected to introduce a significant energy 
dependence on the resulting reaction amplitude. 

We mention here that finally we are left with 38 ad- 
justable parameters in the present work. Their values 
will be given and discussed in the next section. 



V. RESULTS AND DISCUSSION 

The photoproduction reactions of both the neutral and 
charged pious are considered in the present work up to 
a cm. energy of = 1.65 GeV. The free parameters of 
the model as specified in the previous sections are de- 
termined by fitting the available differential cross section 



and photon spin asymmetry data. They are given in Ta- 
bles HirVl Here, we note that in the present work, we 
have attached a momentum cutoff in the loop integral in 
Eq. d?]) of the form A^/(A^-|-g^), where stands for the 
loop momentum and A is an adjustable cutoff parameter 
whose value is also given in Table H] 

The resulting parameters of the effective baryon- 
resonance electromagnetic transition vertices associated 
with the current as explained in the previous section 
are summarized in Table HIl for isospin-1/2 resonances and 
in Table Hill for isospin-3/2 resonances, respectively. We 
emphasize here that the given coupling constant values 
should not be confused with the corresponding physical 
couplings, for they are associated with the current Jj^ 
[cf. Eq. (O] which is only a part of the full current J'' 
given by Eq. (jlSp . More appropriate (physical) resonance 
parameters, i.e., masses and coupling strengths, should 
be associated with the poles of the reaction amplitude, 
and their corresponding residues, in the complex energy 
plane. The pole positions and hadronic residues have al- 
ready been extracted [551] in the Jiilich hadronic reaction 
model. Efforts to extract the resonance electromagnetic 
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coupling strengths from the present model are in progress 
and will be reported elsewhere. 

The other adjustable parameters of the model, i.e., the 
cutoffs in the form factors [Eqs. (jA.17|) and (lA.lSp ] and 
the parameters Ci, C'2, and kq in the nucleon electro- 
magnetic current in Eq. ((23|) are given in Table ITVl while 
the parameters h and (3 in the contact current Af^' [cf. 
Eqs. (HH) and ((301)], are summarized in Table El All ad- 
justable parameters were determined by global fits that 
included all data sets as shown in Figs. \7UM 

The calculated differential cross sections as a function 
of pion scattering angle in the cm. frame are shown in 
Figs. [7HTU1 for 7p — >■ TT+n, — )■ 7r°p, 771 — )■ 7r~p, and 
777. — > Tr'^n, respectively, together with the corresponding 
data. Similarly, the calculated photon spin asymmetries 
are shown in Figs. fTTHTH One sees that the overall agree- 
ment with the recent experimental data is very good. 
Some noticeable discrepancies are seen in both the cross 
sections and beam asymmetries, the latter in — >■ ir^p 
and 7p — >■ 7r+n at higher energies. At this point, it is 
difficult to say whether these discrepancies are due to 
the lack of higher-spin resonances in the present model 



or due to the coupled-channels effects other than those 
from r]N and ttA in the pliotoproduction kernel. In this 
connection, we mention that the KT, channel with isospin 
3/2 has been incorporated into the Jiilich model quite re- 
cently [S^l and the inclusion of the isospin-1/2 KT, and 
KA channels are currently in progress. The inclusion of 
the pN, (jN and strangeness channels as well as higher 
spin resonances such as the Di^ and Fi^ resonances into 
our model requires extra free parameters and will be done 
in future work, which may give us a chance to improve 
the quality of the present description of the experimental 
data, especially, in the ~ 1.6 GeV region. 

Another multi-channel dynamical model available in 
the literature which analyzes the pion photoproduction 
reactions is the EB AC model [s^ . This model, based on 
a unitary transformation method, includes two more res- 
onances than the present model, namely the 1515(1675) 
and Fi5(1680) resonances; it also includes two more chan- 
nels in the intermediate state for the photoproduction 
process, namely the pN and aN channels. On the 
other hand, it analyzes neither the 777 — >■ 7r~p nor the 
777 — 7r*'77 reactions. Although they have considered 
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Figure 12. (Color online) Photon spin asymmetry for 7p — >■ yr'^p as a function of the scattering angle. The first and second 
number in each pair of parentheses correspond to the photon laboratory momentum and the nN cm. energy, respectively. 
Data are taken from Refs. [7i.l75l.l80l]. 



older data than those in the present work, thus mak- 
ing a close comparison difficult, their results are of com- 
parable fit quality to the present model results overall 
for both the 7p — >■ 7r°p and 7p — ^ ir'^n reactions. The 
present results are slightly better at higher energies, es- 
pecially for photon spin asymmetries for the 7p ir'^p 
reaction. Anyway, to achieve the level of the fit qual- 
ity of the EBAC model results, no spin-5/2 resonances 
are required in the present model calculation of the pion 
photoproduction reaction. Further studies are needed to 
understand the role of those resonances. 

At low energies, close to threshold, the pion photo- and 
electroproduction reactions are nowadays completely un- 
derstood thanks to ChPT Any meson-exchange 
dynamical model should, in principle, have built in the 
constraints of ChPT. It is, however, not a simple task 
to account for all the constraints dictated by ChPT and, 
in general, only a few basic constraints are taken into 
account in practice. Indeed, building in the chiral con- 
straints into meson-exchange models is one of the ma- 
jor improvements needed for these models. The Jiilich 
model uses the phenomenological chiral Lagrangian of 



Wess and Zumino [Hsj, supplemented by additional La- 
grangians for the coupling of A, uj, rj, oq, and a [32, 33], 
thus honoring some of the chiral constraints. However, it 
also contains phenomenological form factors which spoil 
these constraints in general. 

Nevertheless, it is interesting to see how the present 
meson-exchange dynamical model performs at low ener- 
gies. In Fig. [T5l we show our differential cross section 
results (solid lines) at low energies close to threshold, to- 
gether with the experimental data. These results, how- 
ever, were obtained by a refit of the parameters h and /3 
of the generalized contact current ([HO)) , leaving all other 
parameters unchanged, since it was not possible to ob- 
tain reasonable fits using the same constants as for the 
global fit, in particular, for the n'^p channel. The refitted 
values of {h, f3) are (12.44, 0.04) for the 7r+n channel, 
(4.13, -1.72) for the Tr°p channel, and (6.58, -0.83) for 
the Tr~p channel, respectively. Comparing with those 
values listed in Table |Vl one sees that these values re- 
fitted to the data close to threshold are quite different. 
This should not be surprising in view of the fact that, in 
principle, these parameters are energy-dependent func- 
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Figure 13. (Color online) Photon spin asymmetry for 7?! — > 7r~p as a function of tiie scattering angle. The first and second 
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Figure 16. (Color online) Differential cross section for 7p — >■ Tr^n, jp — > vr^'p and jn Tv~p as a function of the scattering 
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cm. energy, respectively. The solid curves show the results from the full calculation. The dotted curves are obtained by 
switching off the loop integral [the term proportional to the hadronic amplitude T in Eq. (0]. The dashed curves are obtained 
by switching off the contact current other than the KroU-Ruderman term (the term proportional to the meson charge e^) in 
[cf. Eq. (O]. Data are taken from Refs. [zlS. 



tions that are being treated here as constants for sim- 
plicity. Near threshold, this restriction becomes espe- 
cially noticeable. One sees that the overall agreement is 
reasonable, except for the normalization of the differen- 



tial cross sections for jp — >■ 7r°p ai W ~ 1084 MeV and 
for 7p TT+n at = 1080 MeV. 

One should also keep in mind that the present calcu- 
lation is performed in the isospin basis, with averaged 
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Figure 17. (Color online) Photon spin asymmetry for jp — >■ n'^n, jp — s- Tr^p and 771 — > -K~p as a function of the scattering 
angle. The first and second number in each pair of parentheses correspond to the photon laboratory momentum and the ttN 
cm. energy, respectively. The solid curves show the results from the full calculation. The dotted curves are obtained by 
switching off the loop integral [the term proportional to the hadronic amplitude T in Eq. ((Tjl]. The dashed curves are obtained 
by switching off the contact current other than the KroU-Ruderman term [the term proportional to the meson charge e-rr) in 
(cf. Eq. (I30j]- Data are taken from Refs. [zl, [zl, [s^l • 



nucleon and pion masses of 938.92 and 138.04 MeV, re- 
spectively, which corresponds to a common threshold en- 
ergy of about 1077.0 MeV. This value is about 2 McV 
below the correct threshold value of 1079.1 MeV for the 
7P — >■ TT+n reaction and might explain, at least in part, 
why the present calculation over-predicts the cross sec- 
tion a.tW = 1080 MeV which is just about 1 MeV above 
the correct threshold value. In fact, the dashed curve 
in the corresponding panel represents the result obtained 
with the photo-transition amplitudes FsSJj^ and B'^ [cf. 
Eqs. dl]) and ^] calculated with the averaged mass of the 
nucleon and pion set to the neutron and n'^ mass value, 
respectively, while the hadronic scattering amplitudes are 
still calculated in the isospin basis. One sees that the 
agreement with the data is now improved to some ex- 
tent a.t W — 1080 MeV, but not at a higher energy of 
W = 1112 MeV. Similarly, the common threshold en- 
ergy value is about 3.8 MeV above the correct threshold 
value of 1073.2 MeV for the reaction -fp — > 7r"p. There- 
fore, in this case we expect the calculated cross section 
to under-predict the data. Indeed, as can be seen in 
Fig. [15] (solid curve), this is the case even for the data at 
W = 1084 MeV which is about 10.8 MeV above the cor- 
rect threshold value. The dashed curve here corresponds 
to the result analogous to that for the 7r+n final state 
but with the averaged mass of the nucleon and pion set 
to the proton and tt^ mass, respectively. One sees a non- 



negligible improvement in the agreement with the data. 
In the 771 — t:~p reaction, the common threshold value 
is about 0.8 MeV below the correct value of 1077.8 MeV 
and is much smaller than the corresponding differences in 
the other two reactions just discussed. Therefore, in this 
case we do not expect that such a small energy difference 
will affect much the predicted results, especially, if we 
are at energies not so near the threshold as at 14^ = 1085 
MeV shown in Fig. [131 The dashed curve here is the ana- 
log of those for n'^p and 7r°p for the case of the tt~p final 
state, and the effect is negligible as expected. We thus 
conclude that, apart from the isospin symmetry viola- 
tion effects arising from the mass differences of pions and 
nucleons, the present model works remarkably well even 
at low energies close to threshold. A closer comparison 
with the low-energy data, however, requires a calcula- 
tion in the full particle basis. This is especially true for 
the reaction 7p — )■ n'^p, where it is well known that the 
near-threshold cross section results from cancelations of 
competing mechanisms that yield relatively large contri- 
butions individually and, therefore, making the predic- 
tion of its correct threshold behavior a non-trivial issue 

[Slillii- 

The influence of the loop integral in Eq. (O and the 
generalized contact current other than the Kroll- 
Ruderman term in Eq. ([8]) on the differential cross sec- 
tions and photon spin asymmetries are illustrated in 
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the fit. 



Figs. [16] and [TTl respectively. There, the solid curves 
correspond to the results of the full calculation. The 
dotted curves are obtained by switching off the loop in- 
tegral [the term proportional to the hadronic amplitude 
T in Eq. ([T])]. As one can see, the contribution from the 
loop integral (which provides the effect of the hadronic 
final-state interaction) to the cross sections and beam 
asymmetries is very important showing that it cannot be 
ignored in these reactions. Note that, in this work, we 
have the full T in the loop integral of Eq. ([7]) instead of 
the non-pole amplitude X employed in the earlier feasi- 
bility study of Ref. [49| . The dashed curves are obtained 
by switching off all the terms except the usual Kroll- 
Ruderman contact term in the generalized contact cur- 
rent of Eq. As discussed, these contact terms, 
unique to the present approach, are required to main- 
tain gauge invariance of the reaction amplitude. Leav- 
ing them out, significant effects are seen on the cross 
sections in both the 7p — > 7r+n and 7p — ^ n^p pro- 
cesses; the effect in the — >• 7r~p process is negligible. 
These gauge-invariance-preserving terms also influence 
the beam asymmetries at higher energies in a non-trivial 
way. These interaction-current contributions are seen to 



be just as important here as they were found to be for 
NN bremsstrahlung in Ref. as discussed in the In- 
troduction (see also the Summary, Sec. lVII|) . The results 
for both reactions clearly demonstrate that the gauge in- 
variance required for maintaining the internal consistency 
of the full amplitude in a microscopic model is not purely 
a theoretical issue but is necessary for the description of 
the consistent reaction dynamics (see Ref. ^66] for a more 
detailed discussion on this point). 

Figure [TH] shows the total cross section as a function 
of the t:N cm. energy for the 7p — > ir'^n, jp — > ir'^p 
and jn — )■ 7r~p reactions. As one can see, the agreement 
of the full results (solid curves) with the data — the 
latter taken from Ref. [T^I were not included in the global 
fit — is remarkably good over the entire energy range 
considered. This figure also shows the influence of the 
FSI contributions (dotted curves). Again, they are seen 
to be very important. The influence of the contact terms 
apart from the Kroll-Ruderman term in the generalized 
contact current [Eq. (pO)) ] is also shown as dashed 
curves. Here, a significant effect is seen in the jp — ^ Tr+n 
reaction in the energy range W > 1.3 GeV and in the 
7P — > TT^p reaction around the energy W ^ 1.2 GeV. 
Note for the 7p — ir'^p reaction, the effect on the total 
cross section is largely suppressed as compared to that 
on the differential cross sections (see Fig. [TO]) due to the 
phase space factor sin(6'). 

The effects of the photon couphng to the ttA and rjN 
channels on the cross sections and beam asymmetries are 
illustrated in Figs. [12] and [50] respectively. There, the 
solid curves correspond to the full calculation; the dot- 
ted curves and dashed curves are obtained by switching 
off, respectively, the ttA and rjN channels in the loop 
integral in Eq. ([7]). Note that the threshold energy for 
opening of the ttA channel is at 14^ = 1.37 GeV, while 
for the T]N channel is at = 1.49 GeV. One sees that 
the effect of the rjN channel is practically negligible for 
cross sections, while it shows some influence on the beam 
asymmetries at higher energies. The effect of the ttA 
channel is significant for reproducing both the cross sec- 
tion and beam asymmetry data at higher energies. The 
reason for this is that there is an efficient overlap under 
the loop integral between the generalized contact cur- 
rent and the ttA — ttN transition amplitude in the 
£'31 S31 partial- wave state which peaks around 1.6 
GeV due to the strong couphng of the S'3i(1620) res- 
onance to the ttA channel. There is also a significant 
contribution from the Du — >■ Su partial wave due to the 
511(1535) resonance. Here, it should be noted that, since 
so far the hadronic model parameters in the ttA channel 
have not been constrained by the data, the strong cou- 
pling of these resonances to the ttA channel might be 
an artifact of the model in the hadronic sector. While 
the present photoproduction reaction may help constrain 
some of the 7rA-channel parameters, more conclusive re- 
sults can be expected from investigations of two-pion pro- 
ductions ttN — 7> ttttN and jN — >■ ttttN. 

Motivated by the good agreement of our results with 
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Figure 19. (Color online) Differential cross section for 7p — > yr+n, 7p — )■ 7r"p and 771 Tr^p as a function of the scattering 
angle. The first and second number in each pair of parentheses correspond to the photon laboratory momentum and the ttN 
cm. energy, respectively. The solid curves show the results from the full calculation. The dotted curves and the dashed curves 
are obtained by respectively switching off the ttA loop integral and the rjN loop integral in Eq. ((7]|. Data are taken from 
Refs. ^M^- 



the total and differential cross sections as well as the 
beam asymmetry data, we have extracted the multi- 
pole amplitudes for pion photoproduction and compared 
them with those from the George Washington Univer- 
sity's partial- wave analysis [7]. In Fig. [U the results 
for the multipole amplitudes Mi_|_, E2-, and M2- from 
the present calculation (solid curves) are shown together 
with the results from the SAID analysis The latter 
includes not only the cross sections and beam asymme- 
tries but also the target asymmetries and recoil nucleon 
polarizations as well as some double polarization data 
into its analysis. We see that the agreement between 
the two results for the dominant Mi_|_ amplitude is quite 
good. For the smaller i?2- amplitude, the agreement is 
also reasonable, but for the small M2- amplitude there 
is a considerable disagreement. This illustrates the kind 
of uncertainties one should expect from the present-type 
calculations, even though we reproduce the cross sections 
and beam asymmetries quite nicely. It is clear that in or- 
der to extract more reliable multipoles (apart from the 
dominant ones) from the present model, one needs to in- 
clude more independent observables to further constrain 
the model. Actually, the SAID results are also subject 
to some assumptions in their analysis since at present 
there exists no complete set of data. Indeed, in order 
to uniquely determine the amplitude in the present reac- 
tion, one requires at least eight independent observables 



[S?! (See, also a recent discussion [88|] on this issue). In a 
recent analysis [s^. Workman has also investigated the 
sensitivity of the extracted multipole amplitudes to the 
accuracy of the data used in their extraction. 



VI. UNCERTAINTIES 

Regarding the assessment of theoretical uncertainties 
of our results, this is very difficult to do in a quantita- 
tively reliable manner within the present phenomenologi- 
cal effective Lagrangian approach because of the absence 
of a precise ordering scheme for refining the approxima- 
tions. A procedure was outlined in Ref. [s^l for pion- 
nucleon scattering that allows one to quantify how the 
error margins of the data carry over into uncertainties 
of the extracted parameters of the model approach. A 
similar approach could be used as well to assess the sta- 
tistical errors of the photoproduction reaction. However, 
at present, no quantitative error analysis is available for 
the Jiilich model that we employ here for the hadronic 
final-state interaction and so we must postpone such an 
investigation to future work. 

In addition to the statistical errors, there are system- 
atic uncertainties inherent in all phenomenological effec- 
tive Lagrangian approaches that stem from the imple- 
mentation (or violation) of Lorentz covariance, unitarity. 
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Figure 20. (Color online) Photon spin asymmetry for 7p — >■ n'^n, jp — i- n^p and 771 — > -K~p as a function of the scattering 
angle. The first and second number in each pair of parentheses correspond to the photon laboratory momentum and the ttN 
cm. energy, respectively. The solid curves show the results from the full calculation. The dotted curves and the dashed curves 
are obtained by respectively switching off the ttA loop integral and the 77A'' loop integral in Eq. ((Tjl. Data are taken from 
Refs. 17^,^5,, ,80]. 
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Figure 21. (Color online) Real part and the imaginary part of the multipole amplitudes Mi+, E2- and M2- for — >■ nN 
with isospin / = 3/2 as a function of the nN cm. energy. Scattered symbols are amplitudes taken from the George Washington 
University's partial- wave analysis 



analyticity, and (for photoprocesses) gauge invariance, 
from the truncation of reaction channels and from how 
many intermediate resonances are taken into account. 
For the present approach, we expect the last two sources 
of uncertainties to be most relevant. Plans are under- 
way to address these issues in future applications by in- 
cluding more channels {ujN, KA, KY,) and higher-spin 



resonances [such as Di5(1675) and Fi5(1680)]. However, 
it may well be that a better understanding of the sys- 
tematic errors of phenomenological effective Lagrangian 
approaches can only be obtained by comparing the re- 
sults of different formalisms. As an example, we point 
here to the fact that while the present results, by and 
large, are of a quality similar to that of the EBAC model 
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[39| . we do not need the spin-5/2 resonances employed by 
that model. This is an important finding in itself since it 
places the actual role of these resonances in check. How- 
ever, it may also be an indication that these resonances 
are needed in the EBAC model to make up for some 
basic deficiencies of that model, e.g., lack of Lorentz co- 
variance and gauge invariance. In any case, this points 
to the necessity for further investigating such higher-spin 
resonances to get a better understanding of the corre- 
sponding systematic uncertainties, something we shall do 
in the future. 



VII. SUMMARY 

We have presented results for the neutral and charged 
pion photoproduction reactions within a coupled- 
channels dynamical model in conjunction with the Jiilich 
■kN hadron-exchange model. The photoproduction am- 
plitude in the present approach satisfies the important 
properties of analyticity, unitarity and gauge invariance, 
the latter as dictated by the generalized Ward-Takahashi 
identity. The overall agreement with the cross section 
and beam asymmetry data of these reactions is very 
good in the entire energy range considered. Even at 
very low energies close to threshold, we have shown that, 
apart from the isospin-symmetry violation arising from 
the mass differences of pions and nucleons, the present 
model works quite well, especially, in view of the delicate 
cancelations among various competing mechanisms in the 
7P — > TT^p reaction near threshold. A closer comparison 
with the data close to threshold, however, requires a cal- 
culation in the particle basis. 

The present model includes only the spin- 1/2 and -3/2 
resonances, showing that, within this model, there is no 
obvious indication for the need of higher spin resonances 

— in particular, the Di5(1675) and i^i5(1680) resonances 

— in the energy range of up to 1.65 GeV to describe the 
pion photoproduction cross section and beam asymmetry 
data. In this connection, it is very interesting to extend 
the present calculation to other (spin) observables in pion 
photoproduction. 

The appearance of the terms in addition to the usual 
KroU-Ruderman contact term in the generalized con- 
tact current [Eq. ([5(11) ] is a unique feature of the 
present model. These terms account for the compli- 
cated parts of the interaction current that cannot be 
taken into account explicitly at present. Our results 
show that these terms have significant effects on the 
calculated observables in the present reactions. This 
means it is very import to take into account properly 
the gauge-invariance-preserving interaction current for 
the pion photoproduction processes. The importance of 
this current corroborates similar findings reported for the 
A'^A'^ bremsstrahlung reaction quite recently [6^ , where it 
was found to be crucial in reproducing the KVI data [63| , 
something which had eluded theoretical attempts for a 
very long time. The important point here is that this was 



brought about simply by adding the gauge-invariance- 
preserving current M^, as it is determined here for pion 
photoproduction, as a novel four-point contact-current 
mechanism into the description of A^A^ bremsstrahlung, 
without changing any of the other mechanisms for that 
reaction. This reciprocal consistency between the cur- 
rent mechanisms employed in the two processes clearly 
demonstrates that maintaining gauge invariance of the 
reaction amplitudes, as dictated by the respective gener- 
alized Ward-Takahashi identities of all contributing cur- 
rent mechanisms, is not a purely theoretical issue, but 
an indispensable requirement for a consistent, correct 
description of the reaction dynamics, with direct con- 
sequences for our ability to reproduce the experimental 
data (gsI. |66|. 

Apart from the dominant multipole amplitudes, the 
smaller multipole amplitudes extracted from the present 
model calculation are shown to be subject to considerable 
uncertainties, even though the model reproduces quite 
nicely the recent cross section and beam asymmetry data. 
It is clear that other independent spin observables need to 
be included in the analysis to further constrain the model. 
Note that a unique determination of the multipole am- 
plitudes in pion photoproduction requires, in principle, 
at least eight independent observables [s^l which are not 
available at present. 

As we have emphasized in the previous section, the 
nucleon-resonance electromagnetic transition couplings 
displayed in Tables II-III are not the physical coupling 
values and, as such, they are associated only with the 
present calculation. The appropriate physical electro- 
magnetic couplings should be extracted from the residues 
associated with the poles of the photoproduction am- 
plitude in the complex-energy plane. The work in this 
direction is underway and the results will be reported 
elsewhere. 

Finally, we have considered the t:N cm. energies up 
to 1.65 GeV. This upper limit is set by the limitation 
of the Jiilich hadronic model we employed here for the 
hadronic interactions. To analyze data at higher ener- 
gies, one needs to include higher-spin baryon resonances, 
such as £'15(1675) and i^i5(1680). Also, one needs to per- 
haps also include other meson-baryon channels, such as 
the i^A, ATS, and ujN channels. The ATS channel with 
isospin 3/2 has just been incorporated Q into the Jiihch 
hadronic model and the inclusion of the strangeness chan- 
nels with isospin 1/2 is in progress. The KK and A'E 
channels are expected to play a particularly important 
role in rj photoproduction [16|. 

In summary, the present work provides a comprehen- 
sive treatment of pion photoproduction within a covari- 
ant coupled-channels framework based on phenomenolog- 
ical effective Lagrangians. The details of the approach 
have been constructed [lOl with two main goals in mind, 
namely to preserve the gauge invariance of the current as 
an off-shell condition and to allow for the consistent in- 
corporation of the hadronic final-state interaction. Over- 
all, our results show very good agreement with the data 
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and they, moreover, show that both properties are indis- 
pensable if one wants to provide a quantitatively reliable 
description of the reaction dynamics of pion photopro- 
duction across the entire resonance region. The extrac- 
tion of the electromagnetic transition coupling constants 
for resonances from the residues associated with the poles 
of the reaction amplitude is currently in progress. It is 
also straightforward to extend the present approach to 
the production of other mesons, to strangeness produc- 
tion, and also to the electroproduction of mesons. Tack- 
ling all of these reactions is planned for the near future. 
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Appendix: Lagrangians and form factors 

We list here the Lagrangians and form factors used in 
the present work. 

The hadronic interaction Lagrangians are: 
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Table VI. Coupling constant values of the present model fixed 
from independent sources. The references are shown in the 



third column. 


See text for details, nia-^ — 


1260 MeV. 


gNNn 


13.46 




QNNri 


4. YD 


[311 


QNNp 


3.25 


[71] 


Kp 


6.10 


[71] 


gNNu 


11.76 


[71] 




0.00 


[71] 


SNNai 


gjVJV-n- 


[52] 




■m-rr 


SANn 


2.13 


[90-92] 


ffAJVp 


-39.10 


[90-92] 


g^Tvp 


0.11 


[93-95] 




0.32 


[93-95] 


9~tVP 


0.89 


[93, 94, 96| 




0.25 


[93, 94, 96] 



where A in Eqs. (jA.ip and (|A.2p is the mixing parameter 
of the pseudoscalar (A — 1) and pseudovcctor (A = 0) 
type couplings. In this work, A is taken to be zero which 
means we adopt the pure pseudovector type coupling. 
The coupling constant values in the above Lagrangians 
are given in Table IVII All those values have also been 
used for the hadronic part of the amplitude [13, em- 
ployed in the present work. 

The electromagnetic interaction Lagrangians for the 
nucleon and mesons read 



eiV 
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C^^-n = e[(c)^7r) X TrJgA^, 
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g^y-rrp 



g^yriLj 
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(A.9) 
(A.IO) 
(A.ll) 

(A.12) 
(A.13) 



C^^a, = e F^, [2 (d^TT) X - 2 (^-^Tr) x < 

+ 7VX a^^'Jg , (A.14) 

where e stands for the elementary charge unit, and 
e EE (l + T3)/2 and k ee Kp(l + r3)/2 + k„(1 - r3)/2, with 
the anomalous magnetic moments Hp — 1.793 of the pro- 
ton and Kn — —1.913 of the neutron; F^j^^ = d^A^ — d^A^ 
with A^ denoting the electromagnetic field and a'^" = 
d^ai — d'^al] £ap.Xv is the totally antisymmetric Levi- 
Civita tensor with £^123 _ ^j^g meson-meson elec- 
tromagnetic transition coupling constants in the above 
Lagrangians are given in Table IVII Following Ref. [9^ , 
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they are fixed from the decay [94| of the p and uj meson 
into the ^i:^ and 777 channels, respectively. The signs of 
the coupling constants g^-j^p and g^-nui are consistent with 
those determined from the study of pion photoproduction 
in the 1-GeV region The signs of g^rip and g^r^u: are 
inferred from the flavor SU(3) symmetry considerations 
as used in Ref. [ll]. 

The resonance-nucleon photo-transition Lagrangians 
are 



where p and tob denote the four-momentum and mass 
of the off-shell baryon, respectively. The exponent n is 
taken to be 2 for AA^tt vertex and 1 for NNtt vertex. 
The parameter is determined by fitting to the data 
and it is listed in Table ID 

For the hadronic vertices in the i-channel diagrams, 
the following covariant form factor is included 



„3/2± 
-RNf 



2Mn 

(2) 
9RN-f 



WT'^^^F''''d^N + R.C., (A.16) 



where F^"*"^ = 75 and F^~) = 1; the superscript of Crnj 
denotes the spin and parity of the resonance R. 

For the s-channel NNtt vertex and the u-channel NNtt 
and ANtt vertices, the following covariant form factor is 
employed in our model: 



(A.17) 



A? 



(A.18) 



where a stands for the off-shell meson (a — tt, p, ai); 
q and ma denote the four-momentum and mass of the 
off-shell meson, respectively. The exponent Ua is taken 
to be 1 for a = TT [OO, [HI and 2 for other mesons [92] . We 
use the same cutoff A„ for p, uj and oi mesons in order 
to reduce the number of model parameters. The values 
of the cutoff parameters A^ and A„ are determined by 
fitting to the data; they are listed in Table HI 

Note that the gauge invariance feature of our photo- 
production amplitude is independent of the specific form 
of the form factors. 
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